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Abstract 

We show that the dependence of the charged particle multiplicity on the 
centre-of-mass energy of the collision is, in the String Percolation Model, 
driven by the same power law behavior in both proton-proton and nucleus- 
nucleus collisions. The observed different growths are a result of energy- 
momentum constraints that limit the number of formed strings at low en- 
ergy. Based on the very good description of the existing data, we provide 
predictions for future high energy LHC runs. 



1. Introduction 

Measurements of particle multiplicities constrain the early time properties 
of colliding systems. In the nucleus-nucleus case, these measurements are 
an essential ingredient for the estimation of the initial energy and entropy 
densities and thus the initial conditions from which the system will eventually 
thermalize and the quark gluon plasma be formed. 
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Data collected at RHIC and the LHC for proton-proton (pp) and nucleus- 
nucleus (AA) collisions establish unambiguously that nucleus-nucleus colli- 
sions are not a simple incoherent superposition of collisions of the participat- 
ing nucleons, dnAA/ dr}\ Q > N part - dn pp /dr)\ , and thus that multiple scat- 
tering plays an important role. Further, the possible scaling with the number 
of nucleon-nucleon also does not hold dn^/ dr)\ v=0 <C N coU ■ dripp/ drj\ Q , in- 
dicating that coherence effects among the relevant degrees of freedom at the 
nucleon level are at play during the collision process. 

The collision centre-of-mass energy dependence of the charged particle 
multiplicity in both pp and AA collisions is well reproduced by a power law 
as suggested by models, e.g. the Colour Glass Condensate or the String 
Percolation Model, where coherence effects play an important role. The 
logarithmic growth consistent with pre-LHC data is ruled out. Notably, this 
dependence is stronger [T] in the AA (driAA/dr)\ ~ s ' 15 ) case than for 
pp (dupp/drjl Q ~ s 011 ). Several possible explanations for this difference in 
energy dependence, naively at odds with theoretical expectations, have been 
put forward recently (21 El III El El Ej • 

The String Percolation Model (SPM) [8], the Dual Parton model [9] in- 
cluding parton saturation effects, describes consistently properties of bulk 
multiparticle production pm El El 1131 El El [16] . This framework is closely 
related to descriptions based on the Colour Glass Condensate and Glasma 
flux tubes [IT] for which the same issues have been discussed extensively 

[ibi na 

In this short note, we detail how the SPM can provide for a joint descrip- 
tion of the energy dependence of multiplicity in pp and AA collisions once 
energy-momentum conservation constraints are taken into account. Further, 
we show the resulting dependences on the nuclear species A and number of 
participating nucleon N part to be fully compatible with available data. Fi- 
nally, we give predictions for future higher energy LHC runs and discuss 
generic expectations at high energy. 

2. Mid-rapidity multiplicity in the SPM 

The Glauber model and its generalizations (see [21] for a review) relate 
nucleus-nucleus collisions to collisions of their constituent nucleons. In the 
single scattering limit the average number of participating nucleons (per nu- 
cleon) Na = iVp ai . t /2 behave incoherently and the mid-rapidity multiplicities 
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in nucleus-nucleus and proton-proton collisions are simply related by: 



dn AA 



dy 



Na 



dnPP 



y=0 



dy 



y=0 



This result, which corresponds to the wounded nucleon model 
pected to dominate at sufficiently low centre-of-mass energies, 
energies, multiple scattering becomes important, and 



dn 



AA 



dy 



« 3 - N A ) ■ 



dnPf 



dy 



(1) 

[22], is ex- 
At higher 



(2) 



y=0 



Here, iV^ 3 is the total number of nucleon-nucleon collisions and single scat- 
tering has been subtracted [231 [21] ■ 

Energy-momentum conservation constrains the combinatorial factors of 
the Glauber model at low energy. In the framework of SPM, these constraints 
translate into the sharing of energy-momentum of Na valence strings amongst 
N^ 3 (mostly sea) strings. A possible solution to this problem, the reduction 
of the height of the plateau for sea strings, was pursued in [231 [21]. Here, 
we proceed differently and account for energy-momentum conservation by 
reducing the effective number of sea strings via 



with 



a(y/s) 



1 - 



3\ l + ln(v^7^+l) 



We thus can write 

dn AA 



dy 



N A (N% (V ~ S) - I) 



dnPP 



y=o 



dy 



(3) 



(4) 



(5) 



y=0 



such that the wounded nucleon model eq. rtlj) is recovered for y/s y/so, 
and Glauber result eq. Q follows for y/s ^> \/sq, a{\fs) —> |. 

In the SPM one considers Schwinger strings, which can fuse and perco- 
late [25j [26], 123 123], as the fundamental degrees of freedom. Multiparticle 
production is described in terms of these colour strings which are formed in 
the collision and stretch between partons in the parting nuclei and are thus 
longitudinally extended in rapidity. 
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The multiplicity of produced particles dn/dy is proportional to the aver- 
age number of such strings (twice the number of elementary collisions) N s . 
Thus, for a generic NaNa collision, be it pp or AA, one has 



dn 
dy 



~ *n a ■ ( 6 ) 

y=0 



In the impact parameter plane, the colour content of the strings is con- 
fined within a small transverse area Si = ttTq, with r ~ 0.2 ~ 0.3 fm. The 
strings decay via qq and qq — qq pair production and subsequently hadronize 
to the observed hadrons. In the impact parameter plane, the strings appear 
as disks and with increasing energy-density these disks overlap, fuse and 
percolate leading to a reduction of the overall charge [29] . 

A cluster of n strings behaves as a single string with energy-momentum 
corresponding to the sum of the individual strings and with a higher colour 
field corresponding to the vectorial sum in colour space of the colour fields of 
the strings. In this way, the mean multiplicity (ji n ) and the mean transverse 
momentum squared (pj) of the particles produced by a cluster in the limit of 
random distribution of strings are given by (/x n ) = N s F(rj t )(fi 1 ) and (pj) = 
(pip i) / F{rf) where (fix) and (pxi) are the corresponding quantities for a 
single string. For a random distribution of disks, the colour reduction factor 
F(rf) is 

1 1 — e-v* 

where rf is the impact parameter transverse density of strings (disks) 



r,* = N' = ^-N' t (8) 



with Sn a the area of the impact parameter projected overlap of the interac- 
tion. For a NaNa collision, eq. (IS]) can be recast as 



dn 
dy 



~ F{rf) N s Na , N s Na = N;NF**> • (9) 



The colour reduction factor F(if) results in a slowdown of the increase of 
dn/dy with energy and number of participating nucleons. Note that if nu- 
cleons were to act incoherently, as in the limiting case eq. ([!]), Spt A would be 
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given by the area of a single nucleon S p , while once coherence is accounted 
for, eq. (j2l), Sn a is the overall area of interaction. In general, Sn a depends on 
the impact parameter b of the collision (0 < b < 2Ra) with Na —> 0, Sn a — > 



as b — > Ra and Na —> A, S^ A 
constraints are satisfied by 



S A = rrR 2 A = ixRlA 2 /* as b ->■ 0. These 



S 



Na 
A 



(10) 



with (5 > 0. It follows that 



S NA = nRlA^^' 



and 



„t _ n t N a(V^) Al/3 

Vn a — Vp ly A 71 



A 



/3-l 



(11) 



(12) 



Here, motivated by the scaling limit of the number of vertices in a loop-erased 
random walk [30], we set (3 = 5/3 such that 



Vn a — Vp ly A 



A 



N 



2/3 



(13) 



From eqs. ([I]), (|2]), (|5]), (|9| we can write 



1 dn 
N~ A dy 



dn pp 



y=0 



dy 



y=0 



1 + 



- r 



(14) 



The dependence of the multiplicity on the centre-of-mass collision energy y/s 
is fully specified once the average number of strings in a pp collision N* is 
known. At low energy N*, is approximately equal to 2, growing with energy 
as (y/s/m p ) 2X such that 



2 / /-\ 2A 

r \ /V s 



m r 



(15) 



with i? p the proton radius. 
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3. Data description 



The energy and number of participants dependences implied by eq. ( 14 ) 



can be tested against the available experimental data. Before doing so, two 
remarks on the multiplicity given by eq. (14) are in order: (i) it is for all 
particles, while experimental data accounts only for charged particles; (ii) 
it is for mid-rapidity y ~ 0, while data is collected at mid pseudo-rapidity 
r] ~ and although these coincide, the multiplicity value should be rescaled 
by Jacobean of the y — Y rj transformation. Neglecting the dependence of p?, 
on which the Jacobean depends, on centre-of-mass energy and the (small) 
difference between the Jacobean in pp and AA cases, both rescaling factors 
above can be absorbed into an overall normalization constant k. We thus 
write the charged particle multiplicity at mid-rapidity in the pp case as 



dn pp 
aa ch 



drj 



kFU)N 7 



p ■ 



»7=0 



such that, for the general NaNa case, one has 



1 dn™* 



N A drj 



KF{r j t p )N s v [l + 



r)=0 



i) 



(16) 



(17) 



Eq. (17) depends on three parameters: (i) the normalization k; (ii) the 
threshold scale y/so in a(y/s) in e q. (j4 ); (hi) the power A in eq. (15). We 
performed a global fit to pp data [3111321 1331 1341 1351 1361 137] in the range 
53 < v/s < 7000 GeV, and to AA (AuAu, CuCu and PbPb) at different 
centralities (381 EH] f° r 19.6 < y/s < 2760 GeV. The fitted data sample 



consists of 116 points (19 for pp and 97 for AA). 

The best fit yields the parameter values: k = 0.63 ± 0.01, A = 0.201 ± 
0.003, and ^/sq = 245 ± 29 GeV. One notices that while k and A are tightly 
constrained, the threshold scale ^/so is determined with a sizable associated 
error of ~ 10%. This results from the mild dependence of eq. (17), through 
eq. Q, on ^/sq compounded with the sizable errors in the existing measure- 
ments (see Fig. [2] below). 



Fig. [T] shows a comparison of the evolution of the mid-rapidity multiplicity 
with energy given by eq. (17) with data for pp and AA collisions. 



For reference we quote, in Table [TJ the predicted values relevant for future 
LHC runs. 

The dependence of the multiplicity on the number of participating nucle- 
ons implied by eq. (17) is compared with data in Fig. [2j 
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Figure 1: Multiplicity dependence on y/s. pp data from [3T | l32 l l33l l35 l l36 l 137] 
(circles), CuCu (triangles) and AuAu (stars) from [35], PbPb (star) from [35]. Curves 
obtained from eq. (fTTJ) : (N A = 1, A = 1) for pp (grey line); (N A = 50, A = 63) for CuCu 
(blue line); and (N A = 175, A = 200) for AuAu/PbPb (red line). Color online. 



4. Conclusions 



We have shown that, in the SPM, the power law dependence of the mul- 
tiplicity on the collision energy is the same in pp and AA collisions. The 
slower growth in the AA case is due to finite energy-momentum constraints 
which tamper string creation, and thus the multiplicity, at low energy. In 
the high energy limit, F{rf) — > {rf)~ l l 2 and a(y/s) — > 1/3, eq. (17) can be 
written 



1 dn»* N * 



dr) 



dr> pp 
aa ch 



11=0 



dr] 



ri=0 



Na 
A 



1/2 



N 



1/3 



(18) 



With the obtained fit parameters, this asymptotic result becomes a good 



approximation (within 5% of that given by eq. (17)) for y/s > 500 GeV. This 



indicates that finite energy corrections persist to fairly high energies. 
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PP 


PbPb 


^ (TeV) 


8 


10 


14 


3.2 


3.9 


5.5 


1/Na ■ dn^A/dr, 


7J=0 


6.10±0.14 


6.50±0.16 


7.14±0.18 


8.69±0.14 


9.22±0.15 


10.1±0.2 



Table 1: Predicted multiplicities for pp and PbPb at future LHC energies. 

In this high energy limit, the shape of dnN ^ A as a function of the number 
of participants is energy independent, that is to say that the iVp art and y/s 
dependences factorize. In the SPM, energy-momentum conservation results 
in violations of this factorization and they are the origin of the observed 
discrepancy in multiplicity growth with energy in pp and AA. 

The arguments put forward in this short note can be readily adapted to 
the case of asymmetric (proton-nucleus) collisions. Also in this case, energy- 
momentum constraints are expected to play an important role. 
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